Hybrid Simulations of Reaction-Diffusion Systems in Porous Media

Hybrid or multiphysics algorithms provide an efficient computational tool for combining microand macroscale descriptions of physical phenomena. Their use becomes imperative when microscale descriptions are too computationally expensive to be conducted in the whole domain, while macroscale descriptions fail in a small portion of the computation domain. We present a hybrid algorithm to model a general class of reaction-diffusion processes in granular porous media, which includes mixing-induced mineral precipitation on, or dissolution of, the porous matrix. These processes cannot be accurately described using continuum (Darcy-scale) models. The pore-scale/Darcy-scale hybrid is constructed by coupling solutions of the reaction-diffusion equations (RDE) at the pore scale with continuum Darcy-level solutions of the averaged RDEs. The resulting hybrid formulation is solved numerically by employing a multiresolution meshless discretization based on the smoothed particle hydrodynamics method. This ensures seamless noniterative coupling of the two components of the hybrid model. Computational examples illustrate the accuracy and efficiency of the hybrid algorithm.